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Abstract 

Using the self-dual lattice method, we make a systematic search for modular invariant 
partition functions of the affine algebras of g = A%, A\ + A±, G2, and C2. Unlike 
previous computer searches, this method is necessarily complete. We succeed in finding all 
physical invariants for A<i at levels < 32, for G<i at levels < 31, for C2 at levels < 26, and 
for A\ + A\ at levels fci = < 21. This work thus completes a recent A2 classification 
proof, where the levels k = 3,5, 6, 9, 12, 15, 21 had been left out. We also compute the 
dimension of the (Weyl-folded) commutant for these algebras and levels. 



I. INTRODUCTION 



The subject of this paper is the classification of the partition functions of some con- 
formal field theories. The problem can be stated in general terms as follows. Consider a 
conformal field theory that has an operator algebra decomposable into a pair of commuting 
holomorphic and antiholomorphic chiral algebras, g^ and gn, and a Hilbert space which 
can be written as a finite sum of irreducible representations (Ax,, A#) of cjl x cjr with mul- 
tiplicity N\ L \ R . With each such representation (A^, \r), associate a pair of characters x\ L 
and x\ R ■ Then define the partition functions of the theory as combinations of bi- products 
of characters of the form 



The problem is to find, for a given algebra (gL,9R), all combinations (1.1) such that: 
(PI) Z is invariant under transformations of the modular group; (P2) the coefficients 
N\ L \ R are non-negative integers; and (P3) Nn = 1 if A = 1 denotes the vacuum. 

Any function Z satisfying the modular invariance condition (PI) will be called an 
invariant; if in addition it satisfies the condition that each coefficient N\ L \ R > 0, then it 
is said to be a positive invariant; and finally, if the conditions (PI), (P2) and (P3) are all 
met, then it is considered to be a physical invariant. Clearly any conformal theory must 
have at least these three properties to be physically meaningful. 

Such functions appear in the context of statistical physics, for example in the early 
work of Cardy [1] when he formulated a conformal system on a finite rectangular strip with 
periodic boundary conditions. They also arise in the theory of strings, for example in the 
classical paper by Gepner and Witten [2] where the authors discussed the compactification 
of closed strings on Lie group manifolds. The action for a closed string is the same as that 
for a nonlinear sigma model with properly chosen couplings (i.e. the Wess-Zumino- Witten 
model) and is known to be not only conformally invariant but also invariant under the 
much larger current algebra induced by the underlying group, which forms itself a subclass 
of the Kac- Moody algebras. In the context of the closed string theory the object of interest 
is the factor pertaining to the group manifold of the vacuum-to-vacuum string amplitude 
in the lowest-order (one-loop) approximation. 

It is quite possible that, among all such invariants, some particular ones may lead 
to phenomenologically interesting, even viable, string models, as some recent analyses [3] 
have suggested. Also, it has long been a hope by many that a good understanding of the 
partition functions could make a positive contribution to the broader and as yet incomplete 
task of classifying conformal field theories. This hope is justified partly by the fact that 
a conformal field theory can be identified in some sense through its partition functions, 
and partly by some recent progress in the computing and understanding of the invariants. 
The characters of an algebra, considered in their full dependence on the coordinates of the 
space on which they are defined, are linearly independent [4]. And so are functions of the 
type X\lX\ r - K means that two theories with the same partition functions are identical 
and, conversely, different theories must correspond to different partition functions. 

Many modular invariant partition functions have already been known, and much in- 
sight has been gained through the efforts of many, in particular [5], [6], and [7], among 
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others. But classification proofs, which determine all the physical invariants that belong to 
a certain class, exist only for a few cases, namely, the untwisted Kac-Moody algebra 
[5,8] and the coset models based on it, such as the minimal unitary Virasoro models [5] 
and the N = 1 minimal superconformal models [9]; recently a classification proof for 
has also been obtained [10]. These proofs hold for an arbitrary level (a number through 
which affine representations can be related to finite ones). In addition, the complete list 

of level-one partition partitions for the simple Lie algebras An \ Bn , Cn \ Dn \ and the 
five exceptional algebras has been found [11,12]. 

Systematic numerical searches for invariants of algebras are useful in guiding and 
confirming theoretical analyses, and in discovering new exceptional invariants. Several of 
these searches have already been made. For example, [13,3] use a procedure based on 
Verlinde's formula rephrased as an eigenvalue problem. However in general these searches 
are not necessarily complete - i.e. it is possible that an elusive and presently unknown 
exceptional invariant could escape them (e.g. [13,3] place an upper bound on the sizes of 
the coefficients N\ L \ R ). In the work reported in the present paper, we follow yet another 
method, proposed by Warner [14] and independently by Roberts and Terao [15]. In their 
approach, invariants are generated by using the Weyl-Kac formula of characters and the 
theta functions associated with even self-dual lattices. It has been shown [12] that this 
method is complete: it will yield all the invariants of a given class. It is a number-theoretic 
technique that treats all invariants - exceptional invariants as well as members of infinite 
series - on the same footing. It is quite practical, at least for low rank algebras and small 
levels, as in the cases we considered. 

Its completeness and practicality make this lattice method ideally suited to finish off 
the A*2 classification. In [10] this classification was accomplished for all levels k, except 
for k = 3, 5, 6, 9, 12, 15, 21. These levels survived because the proof in [10] broke down for 
them. In this paper our analysis of A^ includes those remaining levels, thus completing 
the classification for A^\ 

We will implicitly assume throughout this paper that the algebras gL,gn of the holo- 
morphic and antiholomorphic sectors of the theory are identical. This is the situation 
most commonly considered; the problem of asymmetric (heterotic) invariants is addressed 
in [16]. 

In Sec. II we will describe the method. In Sec. Ill we discuss the calculations of the 
invariants for the affinizations of the four rank- two algebras: A 2 at levels < 32; G 2 at levels 
< 31; C2 at levels < 26; and Ai+Ai at levels k\ = ki < 21. The results of our calculations, 
namely the complete set of all physical invariants for those algebras and levels, are given 
in Sec. IV. Our concluding remarks are found in Sec. V. The Appendix describes in more 
detail how we find all of the desired lattices. 



II. THE LATTICE APPROACH 

In this section we will briefly review some basic elements of the affine Kac-Moody 
algebras [17,18] and describe the lattice method of Roberts- Terao- Warner [14,15] that we 
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used in our numerical calculations. We will consider here simple algebras g, but analogous 
statements hold also for g semi-simple. 

Let g be a finite-dimensional simple Lie algebra of rank r, and let aY, • • • , a% be its 
simple coroots. The fundamental weights of g are the vectors j3\ , . . . , j3 r defined by the 
inner products ■ ckJ = Sij . Let p be the sum of all the fundamental weights, p = ^2 Pi- 
The coroots span the coroot space. By coroot lattice of g, denoted by M = M g , we mean 
the integral span of its coroots, a subspace of the coroot space. The dual M* of the 
coroot lattice is called the weight lattice, the set of integral weight vectors. Any vector 
A G M* can be decomposed as A = Ym=\ m ifti where m; G Z. If A is a highest weight, all 
mi, . . . , m r are nonnegative integers. 

Let g = gt 1 ) be the untwisted affine extension of g defined by a central element and 
a derivation. Weight vectors of g may then be denoted A = (A, k, n), where A is a weight 
of g, n an eigenvalue of the derivation, and k an eigenvalue of the central element (called 
the level) . Each unitary highest weight representation of g is associated with an extended 
Dynkin diagram with r + 1 nodes labelled by mo, mi, . . . , m r , where mi, . . . , m r are related 
to a horizontal highest weight vector A in the usual way, and the label of the extra node is 
given by mo = k — X^=i m i a i ' ■> where k is the level of the representation and a( are the 
positive integral colabels of g. The number l + J2 a i * s known as the dual Coxeter number 
and denoted by h y . 

An integrable highest weight representation of g is defined as one with mo > 0. From 
this condition it follows that the level k is also a nonnegative integer since for an irreducible 
highest weight representation of a simple Lie algebra mi, . . . , m r are nonnegative integers. 
It is thus convenient to introduce for each level k = 0, 1, 2, . . . a subset of M* containing 
the horizontal highest weight, A, of each level-/c integrable highest weight representation 
of g: 



i=i 



rrii G Z, < m„ 



*}. 



i=l 



For a given level, the number of horizontal highest weights is finite. In the context of 
a conformally invariant current algebra theory, each A G P+(g,k) is associated with an 
'integrable' primary field. Only integrable fields are physically relevant; all of the other 
possible primary fields are nonintegrable and, as such, irrelevant in the sense that any 
correlation function containing one or more of these fields would vanish identically. We 
will usually find it more convenient to use A + p in place of the horizontal highest weight 

\eP+{g,k). 

The Weyl-Kac formula, which proves a very powerful tool in studying modular invari- 
ance, expresses the character of an irreducible integrable highest weight representation of 
g as a certain sum over the Weyl group of g. Since this Weyl group turns out to be a 
semidirect product of the finite dimensional Weyl group, W(g), and a lattice translation 
group, each summand can be recast in terms of a theta series which takes care of the trans- 
lation group, leaving only the summation over W(g) explicit. The result is the celebrated 
Kac-Peterson formula [4]. 

Generally, given any lattice A, the translate v + A of any vector v G Q <E> A defines a 
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glue class. The theta series of that glue class is given by 

S(v + A)(z\r) d = eM^r x 2 + 27iiz ■ x}. (2.1) 

xEv+A 

Here r G C, and z is a complex vector lying in C <E> A. If the lattice A is chosen to be 
positive definite this series converges, and in fact is analytic for all such z and any r in the 
upper half-plane. Given any Euclidean lattice A and a positive number £, we write A^) for 
the positive definite scaled lattice and A~^ for the corresponding negative definite 
scaled lattice. 

The Kac-Peterson character formula for any A G P+(g, k) can now be written as: 

D 9 (z\r)= £ <w)0{^= + M^)(VWw(z)\t), (2.2b) 

wEW(g) * 

where K = k + /i v is the height of the representation, and e(w) =det(w) G {±1} the 
"parity" of the transformation w G W(g). 

Actually, the character should have been written as X\+ P ( u i z i r ) m its full dependence 
on all three coordinates z G C®M and u, r G C, Im(r) > 0. However, the w-dependence is 
trivial, being an overall exponential, and will not play any role in our study; for this reason 
and without loss in generality we set u = 0. On the other hand, it is important to maintain 
the full ^-dependence; otherwise the characters would cease to be linearly independent, and 
the partition functions of two representations that are complex conjugates of each other 
would become identical. 

So far x 9 \ h h as been defined only for A — p G P+(g, k). However, for later purposes, it 
is convenient to use (2.2a) to define X\ k f° r an V ^ £ M*. The following observation [17] 
permits us to evaluate these X\ ■ F° r any A G M* and k > 0, there are two possibilities: 
either the character vanishes 

x f(z,r)=0 (2.3a) 

for all z and r, or there exists a unique w' G W(g) and unique A' G P+(g,k) such that 
X' + p = w'(X) (mod M^) and 

X 9 > : k (z,T) = e(w')x 9 y k +p (z,r) (2.36) 

for all z and r. We call e(A) the "parity" of A defined as e(A) = if (2.3a) holds, and 
e(A) = e(w') if (2.36) does. In the latter case [A] will denote the unique A' + p = w'(X) 
(modM^ 2 )). 

Let us now recall a few relevant facts about lattices [19,20]. An integral lattice is one 
in which all vectors give integral dot products. An even lattice is an integral lattice in 
which all vectors have even norms. An odd lattice is an integral lattice with at least one 
odd-norm vector. A self-dual lattice A is one which equals its dual A*. A is self-dual iff it 
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is integral and has determinant |A| = 1. One calls a gluing A of some lattice A any lattice 
that can be written as a finite disjoint union of glue classes of Aq. In other words, Aq is a 
sublattice in A of finite index, i.e. such that A/Aq is a finite group. 

In the approach we are following, extensive use is made of Lorentzian, or indefinite, 
lattices. In such a lattice, every vector x may be written (xl] xr), and the inner product 
of any pair of lattice vectors x and x' is given by x ■ x' = xl ■ x' L — xr ■ x' R . In par- 
ticular, let us consider the 2r-dimensional indefinite lattice Mg © Mg~ K ^ that we shall 
call A 9 ' k d = (Mg K ^; Mg K ^), with the height of g, K = k + /i v , being used as the scaling 
factor. Its elements are written as x = (xl; xr), where xl, xr G Mg K \ A gluing A 

of A g > k will look like A = & k [ Xl , . . . , x n ] d = {x + £" =1 %ih I h e Z, x G A 9 ' k }, where 
the vectors xi,...,x n G A 9,k * are called the glues of the lattice A. For example, take 
Xi = (Pi/y/K; f3i/\/K) where, as usual, fa are the fundamental weights of g. Then the 
gluing A D = A»' fc [ ], called the diagonal gluing, can be seen to be both even and 

self-dual. (In the appendix we give a practical test for determining whether or not a given 
gluing is self-dual.) 

For any gluing A of the lattice A 9 ' fc , let us define its level k partition function as 

Z\(g, k)(z L z R \T) d = ^2 t(w L )e(w R ) ^2 expfyrz (jx\ - t*x\)\ (2.4) 

x exp[2iriVK ( w L {z L ) ■ x L - w R (z R )* ■ x R )]/Dg(z L \r)D g (z R \r)* . 

Because A is a gluing of A 9 ' fc , (2.2a) tells us we can write Z^(g, k){zLZ R \T) as a sum of 
terms proportional to 

X 9 x"(z L ,r) X ^(z R ,rr (2.5) 

for \l,\ r G M*. Eqs.(2.3a, b) further tell us that Z\(g, k)(zLZR\r) can be written as a 
linear combination (over Z) of terms of the form (2.5), but with Xl — P and X R — p now 
lying in P+(g, k), that is, 

Z A (g,k)(z L z R \T) d ^ Yl ( N a)x l ,x r X^(^,r) X 3 A ' R fc (^,r)*, (2.6) 

where (Na)x l ,x r are integral scalars. For example, the diagonal gluing Ad has the coef- 
ficient matrix (N Ad ) Xl ,\ r = \\W(g)\\ 5 Xl \ r . (||W(^)|| is the order of the Weyl group of 

<]■) 

If in addition the lattice A is both even and self-dual, its associated partition function 
Z A is modular invariant. This is the key property we are looking for. The RTW lattice 
method suggests we generate all these functions (in the Appendix we describe how we 
did so). The resulting invariants are not necessarily linearly independent, nor are they 
immediately consistent with conditions (P2) or (P3). But it is always possible to deduce 
from them a linear independent set and to construct linear combinations 

Z(z l z r \t) = Y A A Z A (g, k)(z L z R \r) (2.7) 
A 
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that satisfy both (P2) and (P3) (precisely how we do this is discussed at the end of the 
Appendix). Such sums will certainly be physical invariants provided each A included in 
the sum is an even self-dual gluing of K 9,k . The set of the partition functions Z\(g, k) 
corresponding to all even self-dual gluings A of N? ,k is non-empty (containing at least the 
diagonal invariant) and finite (its content being limited by the level and the rank of the 
algebra). As shown in [12] this set is complete: it spans the Weyl-folded commutant of 
(g, k) - i.e. the complex space of all partition functions (1.1) invariant under the modular 
group. Thus the lattice method will allow us to compute all possible, and hence all physical, 
invariants for a given (g, k). 

Incidentally, this RTW lattice method has been generalized in [16] to heterotic modular 
invariants, where it is also proven to be complete. It can also be used to find modular 
invariants for the coset models. 

III. INVARIANTS OF THE RANK-TWO ALGEBRAS 

The above formalism can be applied in principle to any Lie algebra g at any level 
k. However, since the dimension of the corresponding lattices increases with the rank r 
of the algebra and the members of the set P+(g, k) grow rapidly in number with r and 
level k, we will limit ourselves in this work to the four rank- two algebras, A 2 , A\ + A±, 
G 2 , and C 2 , and to levels approximately less than 30. Following the approach outlined 
above, we have calculated the partition functions for all of the even self-dual gluings of 

(i4 (*+3). A ( fc+3)) ^ {A (k + 2) A (k + 2 ); A (k + 2) A ( fc+2)) for k = 1% 

From these results, we deduced the corresponding invariants for G 2 and C 2 . In each case, 
a set of linear independent invariants was obtained and finally, by linear combinations, the 
physical invariants satisfying conditions (P2) and (P3) were determined. 

In the algebra g = A2, the two colabels are equal to 1 , and so the dual Coxeter number 
is h y = 3 and the height K = k + 3 at level k. Let fa and fa stand for the fundamental 
weights of A 2 , and p = fa + fa, their sum. The coroot (= root) lattice is also denoted 
by A 2 , and its dual by A 2 . As usual, any A = mfa + nfa G A 2 is identified through its 
Dynkin labels as (m, n). 

The set of all possible horizontal highest weights corresponding to level k is given by 

P + (A 2 , k) = {mfa + nfa \m,n eZ, < m, n, m + n<k}. (3.1a) 

The trivial representation (the vacuum) at level k defined by the highest weight A = is 
associated with the character Xp = so (P3) reads A"n ii = 1. We will also refer to the 
set of highest weights 

P(A 2 ,k) = P+(A 2 , k) + p = {mfa +nfa \m,n G Z, 0<m,n, m + n < k + 3}. (3.16) 

The first step in our procedure is to find all self-dual gluings of (A 2 ; A 2 K ^), a step 
which we will describe in some detail in the Appendix. Next, to convert a lattice partition 
function (2.4) to the form (1.1), it suffices to write xx f° r arbitrary A G A 2 in terms of x\' 
for A' G P(A 2 , k). This is not difficult once the action of the six elements of W(A 2 ) on A 2 
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is written out explicitly. For example, Xm,n = iff m = (mod K), or n = (mod if), 
or m + n = (mod if). 

Finally, under the outer automorphism /i of the A 2 algebra, the character transforms 

as 

Xh(m,n) = Xn,m- (3.2a) 

If Z is any invariant with coefficient matrix N, we can construct its conjugation, Z c , by 

N XlXr =N XlMXr) . (3.26) 

The conjugation of any physical invariant is also a physical invariant. We can similarly 
define two other invariants, C Z and C Z C , respectively by C N\ L \ R = N h (x L ),\ R and c N^ lXr = 
Nh(\ L ),h(\ R ), but in this algebra one always has C Z = Z c and C Z C = Z. 

The fact that G2 contains A2 as a subalgebra of equal rank gives us a many-to-one 
mapping from the A2 invariants to the G2 invariants (this was also discussed in [14]). In 
particular, let a± (6(2) denote the long (short) simple root of G2; then the colabels are 
a\ = 2 and = 1, and Kac-Peterson formula (2.2a) tells us 

Xm,n — Xm,m+n Xm+n,m? 

(3.3) 

where the character x belongs to G2 at level k, and the characters x belong to A% at 
level k + 1. Given any (k + l)-level A 2 invariant Z in (1.1), we can construct a /c-level G2 
invariant Z by applying (3.3) to Z — Z c — C Z + C Z C . Unfortunately, although this mapping 
preserves (PI), it does not necessarily preserve either (P2) or (P3). Nevertheless, all G2 
invariants can be obtained in this way from A2 invariants. Thus we can read off the G2 
Weyl- folded commutant of level k from the A2 Weyl-folded commutant of level k + 1; the 
physical invariants are then obtained in the usual way. 

Let us now turn to the case g = A\ + A\. The lattice M g here is the orthogonal 
direct sum A 1 © A 1 = Z^ 2 ) © Z^ 2 ). The two fundamental weights of this algebra are 
(3i = (1^/2,0) and (3 2 = (0,1/v^). The vector p is as usual (5\ + fa- The integrable 
highest weight representations of the Kac-Moody algebra g = A\ + A\ have two levels, k\ 
and k 2: one for each A\ factor. Define the set 



P+(Ai + At, ki, k 2 ) = {mfii + n(3 2 \ m,neZ, < m < k u 0<n<k 2 }. (3.4) 

The characters of this algebra, Xm'n 2 ^ are labelled by (m,n). They are just products of 
the characters of A± 

xtt 2 = xtx k n 2 ■ (3-5) 

When k\ = k2 (the case we will be considering), there is a conjugation that inter- 
changes the two A\ factors. It maps a weight (m, n) to (n,m), as in (3.2a). Through this 
mapping, any invariant Z can be associated with three other invariants, Z c , C Z, and C Z C , 
as for A2. However, unlike for A2, these four invariants will all be distinct in general. 

Because C2 contains A\ + A\ as an equal rank subalgebra, the same trick used for G2 
will equally work here. The colabels are a\ = a\ = \. A formula analogous to (3.3) holds. 
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It defines a surjective map from the (k + 1, k + l)-level A\ + Ai Weyl-folded commutant 
to the /c-level C 2 Weyl-folded commutant as before. 

There is a property already discussed in [12] (also found in [21]) that proves to be quite 
useful in practice. Consider any Xl — p, Xr — p G P+(g, k) and A 9 ' k , the Lorentzian root 
lattice of g scaled up in the usual way. We call a positive integer L the order of (A^; A^) 
in A 9 ' k when, for any integer m, (mA^mA^) E A k iff L divides m. For each £ relatively 
prime to L, (£Xl; £Xr) is related nontrivially to (A^; A^), that is, €(£Xl)€(£Xr) ^ and 

N XlXr = e(£X L )e(£X R )N VXLmR] (3.6) 

(e and [A] are defined in (2.3a, &)). See [12] for examples. 

This property has two valuable consequences. All such [^A^]) form a family of 

essentially equivalent representations, so that just one representative of each family need 
be stored. Another practical implication of (3.6) is that if for some £, e(£XL)e(£XR) = — 1, 
then N\ L \ R = for any positive invariant N. 

Call the family ([£\l]; [^A^]) a positive parity family if e(£Ai,)e(£A^) = 1 for all £ 
relatively prime to the order L. By the positive parity commutant we mean the subspace of 
the Weyl-folded commutant consisting of all invariants Z with the property that N\ L \ R ^ 
only for (Xl', ^r) in positive parity families. The positive parity commutant contains all 
positive invariants and so is the only part of the Weyl-folded commutant we need to 
consider. As we shall see, it is generally significantly smaller than the full Weyl-folded 
commutant. 



IV. THE PHYSICAL INVARIANTS 



Most of this section is devoted to listing all the physical invariants we found. Our 
results are summarized in Tables 1 and 2. We give there the dimensions of the Weyl- 
folded and positive parity commutants, along with the numbers of physical invariants for 
each algebra and level we considered. The positive parity commutant was defined at the 
end of the previous section; it necessarily contains all the physical invariants and was the 
subspace we focused on. The tables show the smallness of its dimension P, especially for 
A2, G2 and C2, compared with the dimension D of the Weyl-folded commutant. This 
strongly suggests that a key component of a classification proof for these (and probably 
all) algebras is eq.(3.6). Indeed that is the case with the A 2 classification proof in [10]. 

It seems to be quite difficult to obtain a general formula for the dimension D of the 
Weyl-folded commutant (though much more is known [22] about a related space, the pre- 
Weyl- folded commutant). Tables 1 and 2 suggest the following formulae for D at prime 
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heights K = k + h> 



fnr a . D _((K + 5)/3 when K=l (mod 3 
2 \ (K + l)/3 when K = 2 (mod 3 



for G 2 : D 



(K + 5)/6 when K=l (mod 3 
(K + l)/6 when K ee 2 (mod 3 



far A +A • n-/(^ + 3 )/ 2 when K ee 1 (mod 4 
tOT • ^" \(K + l)/2 when K ee 3 (mod 4 



for C 2 : D = 



(K + 3)/4 when K ee 1 (mod 4 
(K + 1) /4 when K = 3 (mod 4 



Of course all these formulae only apply to prime K; formulae for composite heights will 
be more complicated and are presently unknown (with one exception noted below). The 
above formula for A 2 is proven in [23]. Ph. Ruelle [24] has also obtained the formula for 
G 2 given above. In addition, at the heights K = p 2 of A 2 , where p = 5 (mod 12) is prime, 
E. Thiran [24] has found the formula D = 2(7q 2 - Sq+ 1) where q = (K + 1) /6. For k = 22 
(the only point of intersection of Thiran's formula and our work), both this formula and 
Table 1 gives D = 10. 

One curiosity we obtained was that all invariants, physical or otherwise, for all the 
algebras and levels we considered were symmetric - i. e. their coefficient matrices satisfied 
N\\i = Ny\. This relation is not true of all algebras (e.g. g = D 4n ), but whenever it 
does hold, an easy calculation shows that the coefficient matrices of all invariants will 
commute with each other, so it could be a valuable theoretical property when studying 
the commutant. We have been able to prove, using the analysis in [23], that all invariants 
of A 2 height K will be symmetric whenever if is a product of distinct primes and is not 
divisible by 3 (a similar result will then necessarily hold for these heights of G 2 ), but the 
case of general K seems more difficult. 

Letting P k denote the set P(A 2 ,k) in (3.16) and dropping the superscripts l A 2 ,k' 
from the characters, the known physical invariants of A 2 become: 

Ak= E l^l 2 ; (4-la) 
\eP k 

T>k d = E Xm,nX«fc(™»-«)(m,n)' for k ^ ( m ° d 3 ) and k ^ 4 ? ( 4 ' 16 ) 

(m,n)eP k 

Vk = 3 Yl \x m ,n + Xco(m,n)+X^(m,n)\ 2 for fc ee (mod 3); (4.1c) 

m = n (mod 3) 

£5 = |Xl,l +X3,3| + |Xl,3 + X4,3| + |X3,1 +X3,4| 

+ |X3,2 + Xi,e| 2 + |X4,i + XmI 2 + |X2,3 + X6,i| 2 ; (A.ld) 

ci 1 ) def I 1 1 1 1 , |2 

£9 = 1X1,1 + Xl,10 + Xl0,l + X5,5 + X5,2 + X2, 



.51 



+ 2|X3,3 +X3,6 +X6,3| 2 ; (4.1e) 
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c( 2 ) def I i i 1 2 i I , , |2 , ol |2 

£9 = 1X1,1 +X10.1 +Xl,10| + |X3,3 + X3,6 + X6,3| + 2|x 4 , 4 | 

+ 1X1,4 + X7,l + X4,7| 2 + |X4,1 + Xl,7 + X7,4| 2 + |X5,5 + X5,2 + X2,5| 2 

+ (X2,2 + X2,8 + X8,2)X4,4 + X4,4(X2,2 + X2,8 + Xt^Y, (4- 1 /) 

def 

021 = |Xl,l + X5,5 + X7,7 + Xll,ll + X22,l + Xl,22 

+ Xl4,5 + X5,14 + Xll,2 + X2.ll + Xl0,7 + X7,lo| 2 
+ |Xl6,7 + X7,16 + Xl6,l + Xl,16 + Xll,8 + X8,ll 

+ Xll,5 + X5.ll + X8,5 + X5,8 + X7,l + Xl,7| 2 ; (4.1#) 

together with their conjugations Z c . The subscript in Ak etc. denotes the level of the 
invariant. Note that T>% = "Dg, = Dq, Sg 1 ^ = £^ c , and £21 = £ 2 i- The invariant (4.16) 
was first found in [25], while (4.1c) was found in [26]. The results (4.1d,e,g) were found 
in [27], and (4.1/) was found in [6]. 

This list was proven complete in [10] for all but finitely many levels. The remaining 
levels are all covered by this program, so we now know eqs.(4.1) give the complete list of 
A 2 physical invariants. 

Now letting P k denote the set p + P + (G 2 , k) and again dropping superscripts from 
the characters, the known G 2 physical invariants are the series 

A= £ IxaI 2 , (4.2a) 
\eP k 

two invariants due to conformal embeddings at k = 3 and k = 4 (found by [27]) 

£3= f |xi,i + X2,2| 2 + 2| X i,3| 2 , (4.26) 
4 1] = |Xi,i + Xi, 4 | 2 + |X2,i + Xi, 5 | 2 + 2|X2, 2 | 2 , (4.2c) 

and another invariant at k = 4 (found in [13]) 

c( 2 ) def I 1 2 1 1 1 2 1 1 1 2 1 1 1 2 1 1 |2 
£4 = IXl.ll +IX2,2| +|Xl,3| +|X2,3| +\Xl,4\ 

+ X2.1X1.5 + Xl,5X2,l + X3,lXl,2 + Xl,2X3,l- ( 4 - M ) 

Eqs. (4.2) exhaust all known G 2 physical invariants. This list has been proven complete 
for all prime heights K = k + 4 satisfying K = 5, 7 (mod 12) [24]; of course our program 
verifies its completeness for all other k < 31. 

For the algebra Ai + Ai, we have only investigated the case k\ = k 2 = k = K — 2. 
The physical invariants belonging to infinite series are all simple current invariants [7], 
and have been described in [3]. We shall label them using simple current notation. These 
invariants may be given in terms of their coefficient matrices Nij^rjr. 

In the next few paragraphs, we mean by Ak, T>k, £10, £16, and £ 2 $ the various Ai 
physical invariants [5] . Tensor products of their coefficient matrices are coefficient matrices 
of Ai + Ai physical invariants. These tensor products give rise to some of the series 
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solutions. In all, there are 4 series when k is odd and 12 series when k is even, and a 
number of exceptional invariants at various levels. 

For k odd, the 4 series are: 

(i) the diagonal invariant N k (0) = f A k <S> -4fc; 

(ii) the invariant N k (Ji J2), for > 1, defined by 



Nk{J\J2)ij,i'j> — 



5n'8jj> 

i,K-i'Sj,K-j> 



if i = j (mod 2) 
otherwise 



as well as their conjugations N k (0) c and N k (JiJ2) c 

For k = 2 (mod 4), there are 12 series: 
(i) the diagonal invariant N k (0); 



def, 



(ii) the invariant N k {Ji) = V k Ak for fe > 2; 



(iii) the invariant N k (J 2 ) d = A k ® V k for k > 2; 

(iv) the invariant N k (JiJ 2 ) defined by 

iV fc (JiJ 2 ) i: 



def 



bu'bjji 





+ 8i,K-i'8j,K- 



if z ^ j (mod 2) 
otherwise 



def, 



(v) the invariant N k (Ji; J 2 ) = T> k ®V k for k > 2; 

(vi) the invariant N k (Ji; J1J2), for > 2, defined by 



N k (Jv, JiJ 2 ), 



def 







if z ^ j (mod 2) 



= { Su'Sjj' + 5i,K-i>5j,K-j> if i = j = l (mod 2) ; 
8 it K-i'8jj> + 8u>8j,K-j' if i = j = (mod 2) 



together with their conjugations iV c . We find N 2 (Ji J2Y = N 2 (JiJ 2 ). 

For k = (mod 4), the 12 series are: 
(i) the diagonal invariant iVfc(O); 



def. 



(ii) the invariant N k (Ji) = T> k <g> A 



(iii) the invariant N k (J 2 ) = f A k <8> ffe! 

(iv) the invariant N k (JiJ 2 ) defined by 

def 



N k (JlJ 2 )ij t i 







8a'8jj> + Si } K-i'Sj,K- 



if i ^ j (mod 2) 
otherwise 



def. 



(v) the invariant N k {J\, J 2 ) — V k (g) D^; 

(vi) the invariant N k (aut) defined by 



N k (aut)i 



8u'8jj> 

!' : ! 8i,K-i'8j,j' 

8u'8j,K-j' 

8i,K-i'8j,K-j' 



if i = j = 1 (mod 2) 
if i = 1, j = (mod 2) 
if i = 0, 3 = 1 (mod 2) ' 

if % = j = (mod 2) 



(4.3a) 



(4.36) 



(4.3c) 



(4.3d) 



(4.3e) 
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together with their conjugations N c . 

Besides these series, we also obtain a number of exceptional physical invariants. 

They occur at any level k where there exists an A\ exceptional invariant; invariants of this 
type are given by the tensor products Aio®£\o, £>io<S>£io, £\q®A\q, ^lo^'Cio, £10® £10, 
.4.16 <S> £w, ^16 ® £i6 5 £16 <8> 4-16, £w <8> ^16 an d £16 <S> £16, and their respective conjugations 
Z c . There exist further sporadic exceptional physical invariants not of this type, which we 
denote by £'l to distinguish them from the A\ exceptionals. They are: 



|2 

|2 



Z'l = 1X1,1 + Xi,5 + X5,i + X5, 5 | 2 + 2|xi,3 + X5, 3 | 2 + 2|x3,i + X3, 5 | 2 + 4|x 3 , 3 | 2 ; (4.3/) 

= f 1X1,1 + X3,5 + X5,3 + X7j\ 2 + |X3,3 + X7,l + Xl,7 + X5,5| 2 

+ |X2,4 + X4,2 + Xe,4 + X4,e| 2 ; (4.3#) 

£g = 1X1,1 + Xl,9 + X9,l + X9,9| 2 + |X3,3 + X3,7 + X7,3 + X7,7| 2 

+ (Xl,3 + Xl,7 + X9,3 + X9,7 + X3,l + X3,9 + X7,l + X7,g)X5,5 
+ X5,5(Xl,3 + Xl,7 + X9,3 + X9,7 + X3,l + X3,9 + X7,l + X7,g)* 

+ 2|X5,5| 2 + |X3,5 + X7,5 + X5,3 + X5,7| 2 + |Xl,5 + X9,5 + X5,l + X5,g| 2 ; (4.3/i) 

£"o (1) = f lxi,i + Xi,7 + X5,5 + X5,n + X7,i + X7,7 + Xn,5 + Xn,n| 

+ |Xl,ll + Xl,5 + X7.ll + X7.5 + X5,7 + X5,l + Xll,l + Xlljl 
+ 2|X4,4 +X4.8 + X8.4 + X8,8| 2 ; (4.3l) 

. . ( c\ \ (Jgf i lOi I o 

£lO l ' = 1X1,1 + Xl,7 + Xll,5 + Xll.lll + 1X5,5 + X5.ll + X7.1 + X7,7| 

+ 1X1,5 + Xl.H + Xll.l + Xll,7| 2 + |X5,1 + X5.7 + X7,5 + X7,ll| 2 

+ |X2,4 + X2,8 + X10.4 + Xl0,s| 2 + |X4,4 + X4,8 + X8,4 + X8,s| 2 + 2|X6,4 + X6,s| 2 

+ |X3,1 + X3,7 + X9,5 + X9,H I 2 + 1X3,5 + X3.ll + X9,l + X9,7| 2 ; (4.3j) 

together with the three conjugations £" (2)c , c £" (2) and C £'{ Q ^ C of (4.3j)- 

The invariant (4.3(7) is due to a conformal embedding. The invariant (4.3/) was found 
in [3], while (4.3/i) was found in [13]. The invariant (4.3z) can be understood as due to 
the conformal embedding of A\ + A\ level (10,10) into C 2 + C 2 level (1,1), applied to a 
non-diagonal physical invariant of the latter. The invariant (4.3j) can be written as the 
product of the coefficient matrix of A\q®£\q with N 10 (JiJ 2 ) (see (4.36)). Note that k = 10 
has an incredibly large number of 27 physical invariants. 

Eqs. (4.3) exhaust all the A\ + A\ physical invariants for k\ = k 2 < 21. The only 
other known ones for k\ = k 2 occur at k = 28. 

As before let P k denote the set p + P + (C 2 , k). The C 2 physical invariants include two 
series of solutions, namely, 

A d = f £ |xa| 2 , (4.4a) 
\eP k 

V k d = Yl IXm.nl 2 + Yl Xm,nX*m,K- m -n k > 1 odd, (4.46) 

m odd m even 
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Xm,nXm,n + Xm.nXm^-m-n f° r ^ 6Ven, (4.4c) 

(m,n)eP fe 
m odd 

and four exceptional solutions: 

£3 = f + X3 )2 | 2 + |xi,4 + X3,i| 2 + 2|x 2 , 2 | 2 ; (4.4d) 

47 def I 1 2 

£7 = |Xl,l + Xl,6 + X3,3 + X7,2| 

+ lxi,3 + Xi,8 + X3,4 + X7,i| 2 + 2|x 4 ,2 + X4,4| 2 ; (4.4e) 

^8 = f |Xl,l + Xl,9| 2 + |X3,3 + X3,5| 2 + |X5,5 + X5,l| 2 + |Xl,7 + Xl,3| 2 + |X5,4 + X5,2| 2 

+ |X9,1 1 2 + X9,l(Xl,2 + Xl,s)* + (Xl,2 + Xl,8)Xg,l + IX5,3|2 + X5,3(X3,6 + X3,2)* 
+ (X3,6 + X3,2)X5,3 + IXl,5| 2 + Xl,5(X3,7 + X3,l)* + (X3,7 + X3,l)Xl,5 + |X3,4| 2 
+ X3,4(X7,1 + X7,3)* + (X7,l + X7,3)X3,4 + |X7,2| 2 

+ X7,2(xi,4 + Xi,e)* + (Xi,4 + Xi,6)X7,2! ( 4A f) 

£l2 d = |Xl,l + Xl,13 + X3,4 + X3,8 + 2X5,5 + X7,l + X7,7 + X9,2 + X9,4| 2 - (4.4#) 

Eqs.(4.46, c) are found in [26], (4.4<i, e, g) are conformal embeddings found in [27], and 
(4.4/) is found in [13]. Eqs. (4.4) exhaust all known Ci physical invariants, and is the 
complete list for k < 26. 



V. CONCLUSION 

In this paper we present a numerical search for physical invariants of the rank-two 
algebras A2, A1+A1, G2, and C2 at levels approximately less than 30. Within these limits, 
made necessary by the modest capacity of our personal computer, we have shown first that 
the self-dual lattice method is a useful and practical tool; in particular, it uncovers in these 
algebras and at these levels, all physical invariants, those belonging to the series as well as 
those lying outside the series. Secondly, we have shown that it can fruitfully complement 
a purely analytical approach. For example, it can identify the physical invariants at those 
levels of A 2 which eluded theoretical arguments. It also can uncover regularities and 
suggest approaches useful for these theoretical arguments. Furthermore, not only does it 
determine the dimensions of various commutants at specific levels but it also can suggest 
general formulas applicable to many other, and in some cases all, levels. It is clear that the 
analytical method and the present numerical method, especially when the restrictions on 
the ranks and levels of the algebras are relaxed, can successfully complement each other 
in working towards a solution of the classification of the conformal field theories. 
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Table 1. Dimensions of the commutants and numbers of 
physical invariants for A 2 and G 2 
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G 2 
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N 
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N 
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Z 
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1 
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9 
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/I 


Q 
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O 


Q 

O 
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/I 
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A 
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1U 


e 
o 


7 


Q 
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1U 


ft 
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4 
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1 
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1 
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1 1 
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/I 
4 


4 


4 


1 
1 


1 
1 


1 9 


Q 

y 


4 


4 


ft 
D 


1 
1 


1 
1 


1 ^ 

lO 


1 9 
1Z 


A 
4 


A 

4 


Q 

o 


1 
1 


1 
1 


14 


D 


4 


4 


p: 



1 
1 


1 
1 


1 ^ 


1 1 
1 1 


4 


4 


A 
4 


1 
1 


1 
1 


1 ft 


o 
o 


4 


4 


1U 


1 
1 


1 
1 


1 7 


91 
Zl 


7 
I 


4 


ft 

D 


1 
1 


1 
1 


1 R 


1 9 
1Z 


c; 



A 
4 


ft 

D 


1 
1 


1 
1 


1 Q 


1 
1Z 


/I 
4 


/I 
4 


A 

4 


1 
1 


1 
1 


90 
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o 


/I 
4 


4 


1 9 
1Z 


1 
1 


1 
1 


91 


9R 
ZD 


Q 

y 


c; 



c; 



1 
1 


1 
1 


99 


10 
J. vy 


4 


4 


7 






23 


14 


4 


4 


5 


1 


1 


24 


12 


4 


4 


15 


1 


1 


25 


26 


4 


4 


5 


1 


1 


26 


10 


4 


4 


21 


1 


1 


27 


28 


4 


4 


6 


1 


1 


28 


12 


4 


4 


13 


1 


1 


29 


24 


4 


4 


9 


1 


1 


30 


17 


4 


4 


9 


1 


1 


31 


18 


4 


4 


11 


1 


1 


32 


22 


4 


4 









k = level; D = dimension of the Weyl- folded commutant; 
P = dimension of the positive parity commutant; and iV = number of physical invariants. 



14 



Table 2. Dimensions of the commutants and numbers of 
physical invariants for A\ + A\ and C2 
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O 


A 


O 
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O 
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O 

z 
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Z 


O 

z 
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A 

4 


4 


4 
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z 


V 
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1U 
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lo 
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z 


7 


D 
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4 
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lo 
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lo 
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O 
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O 


Q 
O 


Q 


D 
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4 


/I 
4 


Q 

y 


4 


9 

z 
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07 
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zo 


07 
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4 




z 



z 
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O 


/I 
4 


/I 
4 


D 


9 

z 
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z 


1 9 
±z 


1 7 
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1U 
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Iz 


Q 

y 
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O 
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O 


10 


99 
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4 
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4 


7 
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Z 
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Z 


14 


9n 
zu 
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z 
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z 
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1 n 
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A 
4 


/I 
4 
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1U 


9 
z 
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z 
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ID 




19 
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2 


2 


17 


10 


4 


4 


14 


2 


2 


18 


44 


18 


12 


11 


2 


2 


19 


26 


4 


4 


9 


2 


2 


20 


25 


10 


12 


6 


2 


2 


21 


12 


4 


4 


21 


2 


2 


22 








9 


2 


2 


23 








11 


2 


2 


24 








12 


2 


2 


25 








17 


2 


2 


26 








8 


2 


2 



/c = level; D = dimension of the Weyl- folded commutant; 
P = dimension of the positive parity commutant; and iV = number of physical invariants. 
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APPENDIX 



We describe an algorithm for finding all even self-dual gluings of 0^ = > ^2 ) > 

where is as usual the level. The method can be generalized to other algebras. 

This algorithm involves the gluing construction of lattices [19,28]. One of the powers 
of gluing is the ease with which e.g. self-duality can be verified. A lattice is self-dual iff 
it is both integral and its determinant equals one. A gluing A = Ao[<7i, . . . ,g n ] is integral 
iff its base lattice Ao is integral, its glue generators gi all lie in Aq, and the dot products 
gi ■ gj are all integers. The following fact [28] allows for straightforward computation of 
the determinant |A| of a gluing: 

|A| = |Ao|/||G|| 2 , 

where \\G\\ is the order of the glue group G = A/Ao (i.e. the additive group spanned (mod 
Ao) by the glue generators ^), or equivalently the number of different glue classes x + A 
in A. 

Let K = k + 3 and I = 3K. Let ol\ and j3\ be the first simple root and its corre- 
sponding fundamental weight of A2. By (a, b) we mean the two-dimensional (Euclidean) 
vector (a, b) = a(3i/\^K + bcei/V~K £ A^* and, similarly, (a, b; c,d) the four-dimensional 
(Lorentzian) vector (a(3±/ ' \l~K-\-ba\j ' \[K\ a[3i/ Vk + bai/VK) £ Q* k . It can be shown that 

any vector in equals (a, b) for some integers a, 6; similarly for f2£ and (a, 6; c, d). Any 

even self-dual gluing A of has a maximal left-hand side sublattice A^ defined by the 
set Al = {xl\(xl;0) £ A} and a maximal right-hand side sublattice A^ defined by the 
set A^ = {£_r|(0; xr) £ A}. Both of these lattices are even gluings of and necessarily 
have equal determinants: \A L \ = \A R \. We will construct all A by first constructing all 
possible A L and A R . 

The first step of the algorithm is to generate all even gluings of A^ . Any gluing of 
A^ can be generated by two independent glue vectors (a, b) and (c, d), and, without loss 
of generality, these vectors may be taken to belong respectively to certain sets Si and S2 
to be described as follows. 

The first set, Si, consists of (£, 0) and all the vectors (a, b) such that < a < £, 
< b < K, a divides £, and a 2 /3 + ab + b 2 = (mod K). The second set, S 2 , consists of 
all vectors (0, d) such that < d < K, d divides K, and d 2 = (mod K). 

The possible lattices A^ and A# are then given by all the triples a, 6, d satisfying the 
following conditions: (a, b) £ Si, (0, d) £ S2, b < d, ib is divisible by ad, and (a + 2b)d = 
(mod K). The three congruences (mod K) correspond to the conditions that the glues 
(a, b) and (0, d) have even norms and integral dot products. The remaining conditions 
are designed to force a unique choice of glue generators. There is an obvious one-to-one 
correspondence between a triple a, 6, d and an even gluing of A^ , given by 

i/a K/d 

a, 6, rf-4 K) [(a,6),(0,rf)] d ^ f (J {A K) +i(a, b)+j(0,d)} = \J \J {4 K) +i(a,b)+j(0, d)}. 

i,jez i=i j=i 

It turns out that (a, b) and (0, <i) form a basis for such a lattice. We are actually generating 
more lattices here than we need later on: several lattices may lead to the same lattice sum 
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(2.4) because of the "Weyl-folding" (summing over W(g)) involved in this calculation; 
Weyl-equivalent gluings may be eliminated at this point. 

The next step in the procedure is to choose the base lattice (Ax,; A#) satisfying the de- 
terminantal condition |Ax,| = \^r\- This is equivalent to choosing triples a, 6, d (associated 
with Ax) and a', b', d! (associated with A R ) satisfying ad = a'd' . 

To construct the four-dimensional even self-dual gluings A, we define first the integers 
Li = £/a, L 2 = K/d, L[ = £/a', and L' 2 = K/d' ', and then the vectors 

hi = (2L U -Li) 

h 2 = (-3L 2 - 2Li6/d, 2L 2 + L t b/d ) 
h[ = (2L[,-L[) 

h! 2 = (-3L' 2 - 2L[b'/d', 2L' 2 + L[b'/d' ) 

These satisfy hi ■ (a, b) = h' i - (a', b') = Sn and hi • (0, d) = h\ - (0, d!) = 8 i2 , so h\ and h 2 (h[ 
and h 2 ) form a basis for A^ (A^). Next we find all integers w, x, y, z, where < to, y < a' 
and < x, z < d', such that the two Lorentzian four-dimensional vectors, (hi; wh^ + xh' 2 ) 
and (h 2 ; yh\ + zh' 2 ), have even norms and integral dot products. To each such choice of 
w, x, y, z corresponds the gluing 

A = (Ax; A^ ) [ (hi; wh[ +xh 2 ), (h 2 ; yh[ + zh' 2 ) ] 

= (4 K) ; A ( 2 K) ) [ (00; a', 6'), (00; 0, d' ), (h x ; wh[ + xh' 2 ), (h 2 ; yh\ + zh' 2 ) ] 

Li L' 2 a d 

= U U U U { n k + i(00; a',b')+j(00; 0,d') + l(hi; wh[ + xh' 2 ) + m(h 2 ; yh[+zh' 2 )). 

i=l J = l 1=1 m=\ 

By construction, the resulting lattice is even and also necessarily self-dual: the number of 
different glue classes in it is given by adL[L 2 = 3K 2 . All even self-dual gluings A of 
correspond to some choice of a, 6, c, a', 6', d' and w, x, y, z. 

We will conclude with brief comments on how we find all invariants satisfying (P2) 
and (P3). Using Gaussian elimination on the set {Z\} of lattice partition functions, we 
construct a basis for the Weyl-folded commutant, as well as a basis {Z^\ . . . , Z^} for a 
subspace (the positive parity commutant defined at the end of Sec. Ill) which necessarily 
contains all physical invariants. We are interested in the linear combinations Z = ^ AiZ^\ 
Even when the dimension P of the commutant is large, it is trivial to find the conditions 
on the Ai equivalent to the demand that the coefficients N\ L \ R of Z be integers: e.g. in 
almost all cases we have considered, this condition is simply that each Ai G Z. In our 
notation (P3) becomes the condition that N pp = 1; for our bases, (P3) is always equivalent 
to the choice A\ = 1. Using the remaining positivity conditions N\ L \ R > 0, we then find 
upper and lower bounds for the remaining Ai (the most common ones are < Ai < 1 
or —1 < Ai < 0). We then use the computer to run through the possible values of A^, 
checking for positivity of all coefficients. As we see in the tables, most dimensions P are 
small, but even for large P this final computer search is fast. 
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